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Abstract: We study the anisotropic nature of the Kerr nonlinear response 
in a beta-barium borate (j3-BaB2C<4, BBO) nonlinear crystal. The focus 
is on determining the relevant %^ cubic tensor components that affect 
interaction of type I cascaded second-harmonic generation. Various exper- 
iments in the literature are analyzed and we correct the data from some of 
the experiments for contributions from cascading as well as for updated 
material parameters. We find that the Kerr nonlinear tensor component 
responsible for self-phase modulation in cascading is considerably larger 
than what has been used to date. We evaluate the impact of using such a 
cubic anisotropic response in ultrafast cascading experiments. 
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1. Introduction 

Cascaded second-harmonic generation (SHG) describes the case where the frequency- 
converted second harmonic (SH) is strongly phase mismatched, and this may create a Kerr-like 
nonlinear phase shift on the fundamental wave (FW). Since the early discovery of this cascaded 
self-action nonlinearity cascaded SHG was for a long time overlooked until DeSalvo et 

al. measured a negative nonlinear phase shift from phase-mismatched SHG using the Z-scan 
method [||] (for an early comprehensive review on cascading, see The exciting promise of 
cascading nonlinearities is seen directly from the scaling properties of the Kerr-like nonlinear 
index it generates: «2casc x ~~ d% s /Ak, where c/ e ff is the effective nonlinearity of the quadratic 
interaction. By simply changing the phase mismatch parameter AA: we therefore have a Kerr- 
like nonlinearity that can be tuned in sign and strength. A particularly interesting and attractive 
property is that a self-defocusing Kerr-like nonlinearity n' 7 casc < is accessible by having a 
positive phase mismatch Ak > 0. 

A popular crystal for cascading experiments is beta-barium-borate (j3-BaB20zt, BBO), see 
e.g. [|5|— [1~5|] . It has a decent quadratic nonlinear coefficient, and because the crystal is anisotropic 
it can be birefringence phase-matched for type I (oo —> e) SHG. For femtosecond experiments 
it has the important properties of a low dispersion, a high damage threshold, and a quite low 
Kerr self-focusing nonlinearity; the latter is important for cascading as we discuss later, because 
the material Kerr nonlinearity will compete with the induced cascading nonlinearity. The main 
goal of this paper is to analyze experiments in the literature where the Kerr nonlinear refractive 



index was measured in BBO JT3[-p(i|], and for the fist time extract all four tensor components 
that are relevant to describe the anisotropic nature of the cubic nonlinearity. 

While the anisotropy of the quadratic nonlinearity is extensively used for optimizing the 
phase-matching properties and maximizing the quadratic nonlinear coefficients in nonlinear 
crystals, the case is different when it comes to treating the cubic nonlinearities: the anisotropic 
nature has often been neglected when studying self -phase modulation (SPM) effects from the 
Kerr nonlinear refractive index, in contrast to the case of third-harmonic generation, as evi- 
denced even in early studies Here we show which anisotropic nonlinear susceptibility 
components the experimental literature data represent. The experimental data will be analyzed 
and we correct the reported values for cascading contributions if necessary. We eventually ob- 
tain complete information of all four relevant tensor components for the BBO cubic nonlinear 
susceptibility. Our results show that the tensor component affecting the SPM of the ordinary 

wave, en = Xxxxx = XyyyY' * s ver y we ^ documented, and the corrected literature data agree 
extremely well in the near-IR with the popular two-band model []23|]. The corrected data we 
present shows that c\\ is substantially larger than what has been used so far in simulations in 
the literature, and an important consequence for cascaded SHG is that this reduces the range 
of phase-mismatch values that give a total defocusing nonlinearity. In contrast experiments 
measuring the tensor components that affect the extraordinary wave interaction are more scarce, 
so new measurements are necessary in order to get more accurate and reliable values. We also 
assess the impact of using the anisotropic Kerr response in simulations of cascaded SHG, and it 
turns out that for BBO the optimal interaction angles are such that the anisotropy plays a minor 
role in the wave dynamics. 

As to analyze the experiments properly and rule out any misunderstandings in the definitions, 
the foundation of the paper is a careful formulation of the propagation equations of the FW 
and SH waves under slowly-varying envelope approximation, see appendix These include 
the anisotropic nature of the frequency conversion crystal in formulating the effective cubic 
nonlinear coefficients. With this theoretical background, we then analyze experiments from 
the literature that report Kerr nonlinearities for BBO. When correcting the experimental data 
for cascading contributions it turns out that the most important Kerr tensor component, which 
controls the SPM of the o-polarized FW, is considerably higher than what has been reported 
and used until now. The analysis presented here should help understanding what exactly has 
been measured, and put the results into the context of cascaded quadratic soliton compression. 

2. Cascading background 

In order to understand the importance of knowing accurately the Kerr nonlinear index, let us 
explain first how cascading works. The cascaded Kerr-like nonlinearity can intuitively be un- 
derstood from investigating the cascade of frequency-conversion steps that occur in a strongly 
phase-mismatched medium [Q]: After two coherence lengths 2n/\hk\ back-conversion to the 
FW is complete, and as Ak ^ the SH has a different phase velocity than the FW. Thus, the 
back-converted FW photons has a different phase than the unconverted FW photons, simply be- 
cause in the brief passage when they traveled as converted SH photons the nonzero phase mis- 
match implies that their phase velocities are different. For a strong phase mismatch AkL ^>2n, 
this process repeats many times during the nonlinear interaction length, and in this limit the 
FW therefore effectively experiences a nonlinear phase shift parameterized by a higher-order 
(cubic) nonlinear index casc given as 




Fig. 1. (a) The definition of the crystal coordinate system XYZ relative to the beam prop- 
agation direction k. This convention is in accordance with the IRE/IEEE standard fci]. (b) 
Top view of the optimal crystal cut for type I oo — s> e SHG in BBO, which has (j> = —tz/2 
and 6 = 8 C for perpendicular incidence of an o-polarized FW beam, and the e-polarized 
SH is generated through type I oo — > e SHG. The specific value of the cut angle 6 C depends 
on the wavelength and the desired application. Angle-tuning the crystal in the paper plane 
will change the interaction angle 6. Capital letters XYZ are traditionally used to distinguish 
the crystal coordinate system from the beam coordinate system xyz that has its origin in the 
k- vector propagation direction. 



where a>\ is the FW frequency, n\ and tii the linear refractive indices of the FW and SH, 
Ak = k% — 2k\, and kj = njCOj/c are the wave vectors. The cascaded cubic nonlinearity will 
compete with the material cubic nonlinearity (the Kerr nonlinear index Kerr ), and the total 
nonlinear refractive index change experienced by the FW can essentially be described as 

A " = «2,tot 7 l = ("lease + «2,Kerr) 7 l ( 2 ) 

where I\ is the FW intensity. Due to the competing material Kerr nonlinearity a crucial require- 
ment for a total defocusing nonlinearity n ] 2 tot < is that the phase mismatch is low enough so 
that |«2 casc | > «2 Ken ( tne so-called Kerr limit). When this happens, cascaded SHG can gener- 
ate strong spectral broadening through SPM on femtosecond pulses, which can be compensated 
in a dispersive element. In a pioneering experiment by Wise's group femtosecond pulse com- 
pression was achieved this way [^J. Since the nonlinearity is self-defocusing the pulse energy 
is in theory unlimited as self-focusing effects are avoided [g], and it can even be used to heal 
small-scale and whole-scale self-focusing effects [Qj, Another very attractive feature of the neg- 
ative self-defocusing nonlinearity is that temporal solitons can be excited in presence of normal 
(positive) group-velocity dispersion, which means anywhere in the visible and near-IR. With 
this approach Ashihara et al. [|7|] achieved soliton compression of longer femtosecond pulses. 
Numerous experiments over the past decade have since been motivated by soliton compres- 
sion of energetic pulses to few-cycle duration I I W |l- Since the defocusing limit and the 
soliton interaction depend critically on the n\ Kell . value, knowing an exact value of the n\ Kerr 
coefficient is crucial. 

3. Anisotropic quadratic and cubic nonlinearities in uniaxial crystals 

In a uniaxial crystal, the isotropic base-plane is spanned by the crystal XY axes, and light 
polarized in this plane is ordinary (o-polarized) and has the linear refractive index n . The 
optical axis (crystal Z-axis, also called the c-axis) lies perpendicular to this plane, and light 
polarized in this plane is extraordinary (e-polarized) and has the linear refractive index n e . In a 
negative uniaxial crystal like BBO, n > n e . The propagation vector k in this crystal coordinate 
system has the angle 9 from the Z-axis and the angle relative to the X-axis, cf. Fig. |jja), and 
the e-polarized component will therefore experience the refractive index n e {9) = [n%/ cos 2 9 + 
« 2 /sin 2 0]~'/ 2 , while the o-polarized light always has the same refractive index. 

Degenerate SHG can either be noncritical (type 0) interaction, where the FW and SH fields 
are polarized along the same direction, or critical (type I) interaction, where the FW and SH are 



cross-polarized along arbitrary directions in the (0,<j>) parameter space. In type the FW and 
SH are usually polarized along the crystal axes (9 — or n/2) since this turns out to maximize 
the nonlinearity, and it is called noncritical because the interaction does not depend critically on 
the propagation angles (both nonlinearity and phase matching parameters vary little with angle). 
A great advantage of this interaction is that there is little or no spatial walk off (a consequence 
of the choice 9 = or % /2). In type I 9 and values can often be found where phase matching 
is achieved. This interaction is very angle-sensitive, which is why it is called critical interaction. 
The slowly-varying envelope equations (S VEA) for degenerate SHG are derived in Appendix 
There the quadratic and cubic "effective" anisotropic nonlinear coefficients were symboli- 
cally introduced. Below we show the expressions for these coefficients relevant for SHG in 
BBO, which belongs to the crystal class 3m, and in all cases we consider the reduced numbers 
of tensor components that come from adopting Kleinman symmetry (which assumes disper- 
sionless nonlinear susceptibilities []28|, Ch. 1.5]); thus BBO has 3 independent %^ tensor com- 
ponents and 4 independent tensor components. An excellent overview of the quadratic and 
cubic nonlinear coefficients in other anisotropic nonlinear crystal classes is found in [Q. 

For an arbitrary input (i.e. FW) polarization various SHG processes can come into play (oo — > 
o, oo —¥ e, oe — >• e, oe — >• o, ee — > e, and ee — > o). This includes both type 0, type I and type II 
(nondegenerate SHG, where the FW photons are cross-polarized). Their respective ff e ff-values 
for a crystal in the 3m point group (which also includes lithium niobate) are [ |l9| ] 
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Note that the correct evaluation of the effective nonlinearity requires to take into account the 
spatial walk-off angle p = arctan[tan(0)n 2 /« 2 ] — 9 (for a negative uniaxial crystal) through the 
substitution 9 — > 9 + p in the above equations. Even if p is not insignificant in BBO (around 
2 — 3° for the typical angles used) the difference in c/ e ff is below a few percent, so it is ignored in 
this paper. BBO is usually pumped with o-polarized light as the quadratic nonlinearity is largest 
for this configuration {dn is the largest tensor component, see App. and because the 00 e 
interaction can be phase matched by a suitable angle 9. For such an interaction, the crystal is 
cut so = — n/2 as to optimize the nonlinearity (it is here relevant to mention that dn and ^31 
has opposite signs in BBO). This has the direct consequence that 00 — > o interaction is zero. For 
this reason, when studying BBO pumped with o-polarized light in a crystal with = —n/2, 
Eq. (Q) shows that d e ff = c/31 sin 9 — c/22 cos 0- 

Note that BBO has historically been misplaced in the point group 3, and in addition there has 
been some confusion about the assignment of the crystal axes (where the mirror plane of the 
crystal was taken parallel instead of perpendicular to the crystal X-axis) p9Q. Unfortunately this 
means that even today crystal company web sites operate with d\ 1 as being the largest tensor 
component (in the 3m point group d\\ =0) and supply crystals apparently cut with 0=0 
(because using the point group 3 combined with a nonstandard crystal axes definition means 
that d\\ cos 30 must be maximized, while with the correct point group, 3m, and correct crystal 
axes assignments, ^22 is the largest tensor component and sin 30 = 1 maximizes the effective 
nonlinearity). In order to sort out any confusion, the optimal crystal cut is shown in Fig. |](b), 
and we hereby urge crystal companies to follow the standard notation. 

In Appendix [X] the SVEA propagation equations included also an "effective" third-order 
SPM nonlinearity, %Q{(Oj\ C0j), and cross-phase modulation (XPM) nonlinearity, %jQ {(Of (Ok), 
which due to the anisotropy had to be calculated specifically for a given crystal class and input 



polarization. The SPM and XPM anisotropic cubic nonlinearities for a uniaxial crystal in the 
point group 3m were found in Eqs. (B13), (B14) and (B15) in [pOp. For a type I interaction 
(oo — > e), where the FW is o-polarized and the SH e-polarized, they are 



Z$(ah;ah)=cu (7) 

xj.fi (C02\ COi) = - 4c 10 sin 9 cos 3 9 sin30 + cj 1 cos 4 9 + |c 16 sin 2 29 + C33 sin 4 9 (8) 
Xg S (C0i;ayz) =|cncos 2 0+c 16 sin 2 0+c io sin20sin30 (9) 

Instead for a type ee — > e interaction, as recently considered for cascading quadratic nonlin- 
earities in lithium niobate [|2^] and periodically poled lithium niobate [|3l|,|32[], we have 

ZeffVl^l) =xfs {Gh\B>l) =Zeff Ol^) (10) 

= — 4ciosin 0cos 3 9 sin 30 + c\\ cos 4 9 + |ci6sin 2 20 + C33 sin 4 9 (11) 

(3) 

We have here used the contracted notation c^ m = Xuli' wnere the indices i,j,k,l can take the 



values X,Y, Zand [|33J 

forju: X^l Y -^2 Z ^ 3 

for m : XXX 1 YYY 2 ZZZ 3 YZZ 4 YYZ 5 

XZZ^6 XXZ^l XYY^S XXY —> 9 XYZ ->• (12) 

Note that by using Kleinman symmetry it is assumed that the chromatic dispersion of the non- 

(3) (3) (3) 

linearities is negligible. However, identities like (©1; ©1) = ^ ff (dh; CO2) = Z eff (a>\\0>2) as 
expressed by Eq. ( [1~0| ) have empirically been found not to hold. Specifically, the tensor compo- 
nents they are calculated from turn out to obey slight frequency variations that can be predicted 
by frequency scaling rules, like Miller's rule (see [|34[-|36|] and also later in Sec. |4.4[ ). 

If the Kerr nonlinearity is considered isotropic, we have xje ((Oi',02) = %g^(ffli;ffli)/3 = 
Xtf[(a>z;a>z)/3 [type I 00 ^ e interaction, taking 9 = in Eqs. (R)-(pb], or ^((»i;(0i) = 

(3) (3) x — 1 

XeB (®2> ft) 2) = AJeff ^l) [typ e ee — > e interaction, cf. Eq.dlOh]. These properties underlie 
the parameter B that we used in our previous isotropic model [p7[]. 

— (3) (3) 

A mistake often seen in the early literature on type I interaction is to take Xsvm = #xpm' 

(3) (3) 

this mistake gives a 3 times too large XPM term. As shown above the identity Xspm = Xxpm 
only holds in the type configuration, where the FW and SH have identical polarizations, but 
importantly it is not an identity that is restricted to isotropic nonlinearities as it holds for an 
anisotropic medium as well. Thus, the error made in the past for type I could either come from 
using directly the propagation equations for an isotropic medium and where two pulses with the 
same polarization interact, but it could also come from generalizing type SHG propagation 
equations to type I, and forgetting the XPM properties for cross-polarized interaction. 

4. Measurement of cubic nonlinearities of BBO in the literature 

4.1. Remarks concerning the experimental conditions for the measurements 

Most experiments aiming to measure the cubic nonlinearities have used the Z-scan method [pq], 
which uses a Gaussian laser beam in a tight-focus limiting geometry to measure the Kerr non- 
linear refractive index. The Z-scan method measures the transmittance of a nonlinear medium 
passing through a finite aperture placed in the far field as a function of the sample position (z) 
measured with respect to the focal plane. As the sample is Z-scanned (i.e. translated) through 



the focus of the beam, the lens effect from the Kerr nonlinear index change will change the 
amount of light recorded by the detector; this gives information about the intensity-induced 
nonlinear index change An and thus the n 2 coefficient, defined phenomenologically from the 
general expansion 

An = n I 2 I+n{l 2 + --- (13) 

The intensity is rarely high enough to allow for any contributions but the n 2 term. The super- 
script / underlines that the n 2 parameter is the intensity-dependent nonlinear index, as one can 
also define an electric-field dependent index, typically as An = ln 2 \E\ 2 . In the other sections 
of this paper we use n 2 Kell . to denote the material Kerr nonlinear refractive index and n 2 casc 
the Kerr-like cascading nonlinear index from cascaded SHG; this notation should suffice to 
avoid confusion with the linear refractive index of the SH, ri2- As the result of a closed-aperture 
measurement can be influenced by contributions from multi-photon absorption, the aperture 
can be removed (open aperture scan). The n 2 contributions then vanish, and only multi-photon 
absorption effects remain. Thus it is possible to separate the two contributions. 

There are some issues with the Z-scan method that may affect the measured nonlinear re- 
fractive index. If the repetition rate is too high, there are contributions to the measured n\ from 
thermal effects as well as two-photon excited free carriers [ f39| ] (for more on these issues, see 
e.g. [[£oJ]). Similarly, a long pulse duration can also lead to more contributions to the measured 
«2 Ken than just the electronic response that we aim to model, and in particular the static (DC) 
Raman contribution is measured as well. However, for BBO the fraction of the delayed Raman 
effects is believed to be quite small, so we will assumte that it can be neglected so the measured 
nonlinearities correspond to the electronic response. 

In the following we will analyze measurements in the literature in order to extract the four 
BBO cubic tensor components in Eqs. (^|)-(§). When converting from the nonlinear susceptibil- 
ities Cfim to the intensity nonlinear refractive index we can use Eq. @. There is a small caveat 
here, because the nonlinear index is defined for SPM and XPM terms, where only two waves 
interact through their intensities, while the c^,,, coefficients are generally defined for four-wave 
mixing. Thus, it does not always make sense to define a nonlinear refractive index (consider 
e.g. the cio component) until the total effective susceptibility is calculated, e.g. through Eqs. 
(0)-(|ll]). Thus, it is safer to keep the susceptibility notation as long as possible when calculat- 
ing the effective nonlinearities. The exception is when measuring the nonlinear refractive index 
of the oooo SPM interaction because the relevant nonlinear tensor component that is measured 
is cn, see Eq. (Q), and thus the connection between c\\ (i.e. Xooio) and n 2 Ken . is unambiguous. 

In the experiments in the literature the cascaded quadratic contributions were often forgot- 
ten or assumed negligible. The exact relations for the evaluating cascading nonlinearities are 
derived in Appendix |b| We will address these issues below by in each case calculate where 
possible the cascading nonlinearity. 

4.2. The c\ \ tensor component 

In this section we present measurements using o-polarized input light, which means that 
the en tensor component is accessed by measuring the SPM component of the o-polarized 

13) 13) 

light %lJ{(0\\(Oi) = Xoooo(— (0\,(0\,— (0\,(0\). The connection to the Kerr nonlinear index is 

w 2,Ken ( W i; W = 3 *eff ( W l > W l )/ 4 »l £ 0C = 3c u /4n 2 e C. 

To our knowledge, the first measurement of any of the cubic susceptibilities was performed 



by Tan et al. [J13J]. They investigated the cascaded nonlinearity in a 10 mm BBO crystal cut at 
6 = 22.8° for phase matching at 1064 nm (they claim also that the crystal is cut with <p = 0, 
which is probably due to the confusion related to placing BBO in the point group 3 instead of 



3m, cf. the discussion on p. ||). Pumping with 30 ps o-polarized pulses from a Nd: YAG mode- 
locked laser, they tuned the BBO crystal around Ak = in a type I oo — > e SHG interaction, 
and by recording the impinging and generated pulse energies vs. angle they found by a fitting 
analysis the coefficient g = 220, where g = ^ /d£ s was the only free parameter in the fit. 
Using c/ e ff = 2.04 pm/V the intrinsic (material) nonlinear susceptibility can then be found as 

Xtat = 9 - 2 x 10 22 m2 / v2 - This ^eff-value was calculated using 9 = 22.8°, = -n/2 and the 
nonlinear coefficients in App. R at 1064 nm. They also performed a Z-scan measurement, and 

- (3) - (31 

the setup was calibrated towards a BK7 glass sample. They found j int / Xbki — ^ -4- ^ v usm g the 
value they propose Xkyli = 4.5 x 10~ 22 m 2 /V 2 , one therefore finds xj® = 6.4 x 10~ 22 m 2 /V 2 . 
We report these results with a bar because their definition of the cubic susceptibility is different 
than what we use: they define the SPM polarization as = £ o^im' ^\^\ 2 /2, which replaces 
our definition Eq. (p3|), so the connection between them is %Q = 3/2x^ (dhl COi)- To check 
this, from the literature we find for BK7 the value n 2 Ken . = 3.75 ± 0.3 • 10- 20 m 2 /W measured 
at 804 nm by Nibbering et al. [^TJ, which corresponds to % (3) = 3.03 x 10~ 22 m 2 /V 2 . This 
is precisely a factor 3/2 smaller than the value Tan et al. mentions, and this confirms the re- 
lationship between the cubic susceptibilities. We now apply Miller's rule to the Nibbering et 
al. result to get for BK7 %^ = 2.90 x 10~ 22 m 2 /V 2 at 1064 nm, so we get a corrected value 
xjnt = 6.08 x 10~ 22 m 2 /V 2 . In our notation their cascading measurements gave 

c u = 6.12±0.53-10~ 22 m 2 /V 2 , A = 1.064 /im (14) 
corresponding to Kerr ((Oi ; (Oi ) = 6.32 ± 0.55 • 10 -20 m 2 /W. Their Z-scan measurement gives 

c u = 4.05 ± 0.52 • 10" 22 m 2 /V 2 , X = 1.064 ^m (15) 
corresponding to «2.Kerr 

{coucoi) = 4.18 ± 0.54 • 10" 20 m 2 /W. The error bars are indicative as 
they do not report the errors in their measurements, but we assume 5% error on g and conser- 
vatively estimate an error of 10% on the Z-scan result. 

Hache et al. @] used a 1 mm BBO cut at B c = 29.2° with 800 nm 100 fs pulses from a MHz- 
repetition rate Tksapphire oscillator. The pump was o-polarized allowing for oo — > e cascading 
interaction. They carried out Z-scan measurements both close to the phase-matching point 8 C 
as well as "far from SHG phase matching". In the latter case they found 

"2,Kerr(«i;«i)= 4 -5±l-0-10- 20 m 2 /W, A = 0.8 jUm (16) 

which corresponds to cu =4.36 ±0.97- 10~ 22 m 2 /V 2 . They claim that cascading contributions 
did not contribute for this measurement, but whether this is true cannot be judged with the 
information at hand, as they did not specify the angle they used for this measurement. 

DeSalvo et al. [ |T6| | used the Z-scan method with a single-shot 30 ps pulse at 1064 nm, which 
had k parallel to the optical axis c, so 6 = 0. Thus, (a) the pump is always o-polarized no 
matter how the input beam is polarized, and (b) the input polarization determines the angle 
(p. Since the specific orientation of the crystal with respect to the input polarization was not 
reported, we take it as unknown. The relevant Kerr contribution at the pump wavelength is 
oooo interaction, which is given by Eq. (^) as simply one tensor component c\ \ . The experiment 
reported «f = 1 1 ± 2 ■ 10" 14 esu at 1064 nm, where the nonlinear index change is defined as 
An = \n\ |£| 2 , and using Eq. (CI) in ^ this corresponds to =2.70± 0.49- 10~ 22 m 2 /V 2 . 
The cascading contributions are from oo — >• e and oo —> o processes. Since 0=0 they have the 
same phase mismatch values as e-polarized light in this case has the same refractive index as 
o-polarized light. At 1064 nm the value is Ak 00 e — kk 00 o = 234 mm -1 . Taking into account 



the two cascading channels Eqs. (Q) and ([|), and using Eq. (j5l|), we get the total contribution 
Xcll = -[sin 2 (30)+cos 2 (30)]1.95- 1(T 22 m 2 /V 2 = 1.95 • 10 22 m 2 /V 2 , i.e. independent 
on the propagation angle (j). For these calculations we used dn = —2.2 pm/V. The d e g values 
are typically reported with a 5% uncertainty [^2p, giving a 7% uncertainty on the cascading 
estimate. This means that when correcting for the cascading contributions we get 

en =4.65 ±0.51 • 10~ 22 m 2 /V 2 , A = 1.064 (17) 

corresponding to r^^CO^COi) = 4.80 ±0.53 • 10~ 20 m 2 /W. At 532 nm nf = 21 ±4- 10" 14 
esu was measured, corresponding to x' 3 ' = 5.22 ±0.99 • 10~ 20 m 2 /V 2 . The cascading is smaller 
here, using Ak ooe = Ak aoo = 1,933 mm -1 we get XcaL = —0.61 • 10~ 22 m 2 /V 2 (where we used 
dxi = 2.6 pm/V as measured at 532 nm [^|), so 

en = 5.82±0.99- 10~ 22 m 2 /V 2 , A =0.532 jUm (18) 

corresponding to n[ Kerr (ct>i ; COj ) = 5.87 ± 1.00- lO" 20 m 2 /W. At 355 nm «f = 14± 3 ■ 10" 14 
esu was measured, corresponding to x^ = 3.54 ±0.67 • 10~ 22 m 2 /V 2 . The cascading is small, 
using Ak ooe = Ak aoo = 11,736 mm -1 we get Xeub = -0.38 • 10~ 22 m 2 /V 2 (where Miller's 
scaling was used to get the value ^22 = —4.4 pm/V), so 

en = 3.92 ± 0.68 -10" 22 m 2 /V 2 , A =0.355 jUm (19) 

corresponding to n' 2 ^ en ((Oi', ft>i) = 3.81 ± 0.74 • 10~ 20 m 2 /W. We must mention that the cas- 
cading value is quite uncertain: the SH wavelength is here 177.5 nm, which is very close to the 
UV poles in the Sellmeier equations (the fit in [^3[| puts the poles at 135 nm and 129 nm for o- 
and e-polarized light, respectively). Finally, at 266 nm they measured nf = 1 ±0.3 • 10~ 14 esu. 
Here cascading is estimated to be insignificant (certainly we cannot estimate accurately it using 
the Sellmeier equations as the SH lies right at the UV pole), so converting to SI we have 

en = 0.26 ± 0.078 • 10~ 22 m 2 /V 2 , A = 0.266 jum (20) 

corresponding to n' 2 KeTI ((0i ; ©1) = 0.24 ± 0.07 • 10~ 20 m 2 /W. 

In a similar experiment, Li et al. [ |l8| ] pumped a 5 mm Z-cut BBO crystal with 150 fs 780 
nm pulses from a Tksapphire oscillator. The propagation was along the optical axis, implying 
= and thus that for both cascading and Kerr nonlinearities the angle <p should not matter as 
the pump always will be o-polarized. Nonetheless, they measured using the Z-scan method two 
different results when flipping the input polarization, namely «j = 4.0 ±0.5 • 10~ 20 m 2 /W when 
polarizing light along the [1 0] direction (X-axis, i.e. <j> = 0), and n' 2 = 3.2±0.5 • 10~ 20 m 2 /W 
when polarizing light along the [0 1 0] direction (Y-axis, i.e. <j> = 7t/2). Cascading was ignored 
as it was considered too small, so let us asses whether this is true. The cascading contributions 
are from 00 — > e and 00 —> o, both having the same Ak ooe = Ak 000 = 55 2 mm -1 . However, 
only one come into play in each case: for = the only nonzero nonlinearity is d'^° = — c/22 
while for (j) = n/2 the only nonzero nonlinearity is d'^ e = —^22- Thus, they turn out to have 
the same cascading contribution, which from Eq. ([[]) becomes «2casc = ~ 1 -30 • 10~ 20 m 2 /W. 
When correcting the measured nonlinearities with this value we get 

n 2,K e ir( fi, i; < °i) =5.30±0.51 ■ 10~ 20 m 2 /W, [10 0], A =0.78 jum (21) 
n 2,Keir(<»i;<»i) =4.50±0.51 ■ 10~ 20 m 2 /W, [010], A =0.78 /im (22) 

corresponding to c n = 5.18± 0.50 • 10~ 22 m 2 /V 2 and c n = 4.39± 0.50 • 10~ 22 m 2 /V 2 , respec- 
tively. When reducing the repetition rate from 76 MHz to 760 kHz they saw similar results. 



Ganeev et al. [ |20| ] used a 2 Hz 55 ps 1064 nm pump pulse and a BBO crystal with = 51°, 
The procedure they used was to calculate the "critical" Ak c where defocusing cascading exactly 
balances the intrinsic focusing Kerr, i.e. where |«2casJ = n 2Kerr By angle tuning well away 
from this point (so Ak ^> Ak c ) they tried to make cascading insignificant. We estimate this to be 
true: specifically, the pump was o-polarized and the crystal was cut with = — n/2 to optimize 
SHG [j44[ ]. Then the only cascading channel is oo — s- e, with Ak ooe = —655 mm -1 , and using 
Eq. (|lj) we get «2casc = 0-^1 ' 10 20 rn 2 /W that is therefore self-focusing (as they also note 
themselves). Using the Z-scan method they measured n 2 = 7.4 ±2.2 • 10~ 20 m 2 /W with the 
uncertainty estimated to be 30%, so correcting for cascading 

»2.Kerr(<»i;<»i) = 7.08 ± 2.22 • 10~ 20 m 2 /W, A = 1.064 jum (23) 

corresponding to c\\ = 6.87 ± 2.15 • 10 -22 m 2 /V 2 . At 532 nm they measured n 2 = 8.0 ± 2.4 ■ 
10~ 20 m 2 /W. Here cascading is more significant as Ak ooe — —194 mm -1 , c/ e ff = 1-67 pm/V 
and «2 casc = 2.76 • 10~ 20 m 2 /W. Therefore the Kerr value corrected for cascading becomes 

»2,KerrOi;<»l) = 5.24 ±2.41 • 10" 20 m 2 /W, A = 0.532 fim (24) 

which corresponds to en = 5.20 ±2.39 • 10~ 22 m 2 /V 2 . 

Moses et al. fll"5] ] pumped a BBO crystal at 800 nm with 110 fs o-polarized light from a 
TkSapphire regenerative amplifier. By angle-tuning the crystal they achieved zero nonlinear 
refraction: this happens when the cascading from the oo — > e interaction exactly balances the 
Kerr nonlinear refraction, which was found to occur at A£ = 3 1 ± 5 n/ mm for a wide range of 
intensities. Using this phase mismatch value and that n\ tot = then a reverse calculation though 
Eq. ([!]) gives n 2 Kcrr = — n 2 casc . Using the value c4ff =1-8 pm/V they inferred the value n 2 Kerr = 
4.6 ± 0.9 • 10" 20 m 2 /W (the main sources of the uncertainty are determining the precise phase- 
mismatch value as well as uncertainty on the c/ e ff used). Since o-polarized pump light was used, 
this contains only the c\\ tensor contribution from an oooo interaction. However, the value 



deft =1-8 pm/V they used to infer this nonlinearity was probably taken too low |45|. Let us 
therefore estimate it again using a more accurate value: Ak = 31 7r/mm is achieved at 800 nm 
with = 26.0°, giving d e s = 2.11 pm/V (found using a?22 = —2.2 pm/V and ^31 = 0.04 pm/V 
at 1064 nm and employing Miller's rule to convert to 800 nm). With this corrected value we get 

»2,Kerr(®i;®i) =5.91 ±1.04- 10" 20 m 2 /W, A =0.8 jum (25) 

which corresponds to en = 5.77 ± 1.02 • 10~ 22 m 2 /V 2 . 

This completes the en measurements. A summary of the data is shown in Fig. g, together 
with the predicted electronic nonlinearity from the two-band model (2BM) p3[|. The experi- 
mental near-IR data are amazingly well predicted on an absolute scale by the 2BM; we must 
here clarify that for the 2BM we chose the material constant K = 3100 eV 3 / 2 cm/GW, which 
was found appropriate as a single material parameter for dielectrics [|6j], and by using the BBO 
band-gap value E g = 6.2 eV []16|]. Thus, in practice there are no free parameters in the model, 
which underlines the incredible agreement obtained on an absolute scale (and not just the trans- 
lation from one wavelength to another). We mention here that in DeSalvo et al. [|l6|| they pro- 
posed to "rescale" the bandgap do obtain a better fit with the experimental data (mainly because 
the two-photon absorption values /3 were not accurately predicted with 6.2 eV), and eventually 
suggested using 6.8 eV instead of 6.2 eV for BBO. This gave a much better j3 agreement at 
266 nm, while the 352 nm value was still off. Here we see that the corrected DeSalvo n 2 j^ err 
near-IR values actually agree extremely well with the 2BM nl, Kerr value when using E g = 6.2 
eV, while the agreement is less accurate with E g = 6.8 eV. A common issue whether one uses 




Fig. 2. Summary of the experimental data for n ! 2 -values from the literature corresponding 
to the en nonlinear susceptibility coefficient (n' 2 = 3c\ \ /4n 2 £pc). The plotted values are 
the ones reported in Sec. 4.2. References: Tan et al. 1993: B13|]; Hache et al. 1995: [14]; 
DeSalvo et al. 1996 Li et al. 2001 Ganeev et al. 2003~Q]; Moses et al. 2007 [15]. 
The theoretically predicted electronic nonlinearity is calculated with the 2-band model [ 23 ] . 
The average value curve was calculated through a weighted mean of the Miller's delta from 
all data, except the UV measurements below 400 nm, and the shaded areas represent one 
and two standard deviations. 



6.2 or 6.8 eV is that the predicted strong enhancement just below 400 nm due to two-photon 
absorption of the FW is not reflected in the experimental data. This is of little importance for 
our purpose, because in cascaded SHG the FW wavelength is usually never below 600 nm. In 
order to extract a suitable value for the cn coefficient, we calculated first the Miller's delta from 
all the data. It is defined as 

(3) 

xPx^x^x^ 

where %, (do) = w y(do) — 1 is the linear susceptibility. We here excluded the UV measure- 
ments below 400 nm as they obviously are not represented well by Miller's scaling. Since 
Miller's delta should be independent of wavelength this should allow us to calculate an average 
value based on the measurements at different wavelengths. We confirmed this hypothesis by 
checking that that Miller's delta vs. measurement wavelength could be fitted to a line with a 
near-zero slope, before we calculated the weighted average (with a weight given by oj 2 where 
Oj is the estimated error of the y'th measurement) as 

A U11 =52.8 ±8.4 x 10~ 24 m 2 /V 2 (27) 

where the standard deviation was calculated using an unbiased weighted average. This value 
corresponds to c u = 5.03 ■ 10~ 22 m 2 /V 2 and n^ Kerr (di;di) =5.15 ■ 10~ 20 m 2 /W at 800 nm. 
We see in Fig. |^ that this average is a quite good representative of the experimental data; all 
except one are within two standard deviations. 



4.3. Other tensor components 



Banks et al. [19] measured the cio tensor component of BBO using third-harmonic generation 
(THG) measurements at 1053 nm using 350 fs pulses from a 10 Hz TkSapphire regenerative 
amplifier. The propagation angle was 9 = 37.7°, and from the data recorded while was varied 
around zero they extracted Cio = —6 • 10 -24 m 2 /V 2 using the quadratic nonlinear coefficients 
of Shoji et al. [|42"|]. (They extracted other values as well, but we choose to use this one because 
it was calculated from the experimental data using the same effective quadratic nonlinearities 
of Shoji et al. as we use.) In fact, Banks et al. reported that this value came with a quite large 
uncertainty due to the uncertainty in the d$\ - and dis, in absence of Kleinman symmetry - 
coefficients in the literature. Since [[1^] define C^,„ = %^ /4, we have C^ m = c^ m /4, so 

ci = -0.24 ± 0.04 -10" 22 m 2 /V 2 , A = 1.053 Aim (28) 

The uncertainty of the measurement was not reported, but for a similar fit it was reported to 
be 15%, which is what we used above. They also measured a three times larger value by using 
different quadratic nonlinearities, namely d^ = 0.16 pm/V (as opposed to d\s = 0.03 pm/V). 



Banks et al. [|19| also measured the mixture iCn cos 2 9 m + C^sin 2 6 m = 4.0 ± 0.2 



10~ 23 m 2 /V 2 at 9,„ = 47.7°. If we use the Moses et al. result, appropriate ly s caled to 
1053 nm and for CO + CO + co — > 30) interaction using Miller's scaling (see Sec. |4.4| ) to give 
en = 6.07 • 10~ 22 m 2 /V 2 , then we get 

c 16 = -1.23 ±0.36 -10" 22 m 2 /V 2 , A = 1.053 fim (29) 

A small caveat must be noted: what they measured was related to the %^ (— 30); CO, CO, co) coef- 
ficients, because they investigated the yield of the third harmonic with respect to the pump, 
and decoupled the cascade yield (from multistep SHG mixing, i.e. co + CO — > 2co followed 
by CO + 2co — > 30)) to get the pure cubic nonlinear contribution. However, it is not sure that 
X^\—3co;co,co,co) = %^\—co;co,-co,co), i.e. that the THG nonlinearity is the same as the 
SPM nonlinearity; this was discussed, e.g., in [^o|. Therefore using these values to model the 
SPM and XPM effects in BBO is an approximation. 

Sheik-Bahae and Ebrahimzadeh [ |l7| ] used the Z-scan method to measure n' 2 = 3.65 ± 0.6 • 
10~ 20 m 2 /W (£ (3) = 3.43 ±0.56 ■ 10~ 22 m 2 /V 2 ) using a 76 MHz TkSapphire oscillator giving 
120 fs pulses at X = 0.850 jim. The BBO crystal was cut for phase matching at 800 nm, i.e. 
with 9 = 29.2°. We assume the crystal cut is optimized for SHG, so <p = —n/2 (which makes 
sin30 = 1). The pump pulse was e-polarized as to ensure lack of phase matching and thereby no 
cascading. We will now assess whether this assumption is fulfilled. The cascading contributions 
that might occur are ee — > o and ee — > e, however the former has zero nonlinearity when = 
—n/2, cf. Eq. (|5|). The phase mismatch for the latter interaction is Ak eee = 396 mm -1 for 9 = 

29.2°. The cascading contribution at = —n/2 then becomes Xcllc = —0.69 • 10~ 22 m 2 /V 2 . 
Thus, the choice of using e-polarized light does make the cascading contribution small, but not 
insignificant (it is roughly equal to the reported uncertainty). If we now correct the measured 
value with the cascading contributions we get Xk,L = 4-08 ± 0.57 ■ 10~ 22 m 2 /V 2 . Let us now 
understand what tensor components this value represents. By using e-polarized light several 
tensor components come into play for the cubic nonlinearity: the effective Kerr nonlinearity 
experienced by the pump is an eeee interaction, which is given by Eq. (JIT]). For 9 = 29.2° and 
(j) = —n/2 the effective nonlinearity is 

X$ (»i ; coi ) = - 1 .30cio + 0.58ci i + 1 .09ci 6 + 0.057c 33 (30) 

Four tensor components appear, and we note that we here have a measurement that involves 
the C33 component. Let us now use the previous results to extract it: using Eq. ([251), (jH) an d 



Table 1. Summary of the literature measurements of the cubic self-action (SPM or XPM) 
nonlinearities. The column (A) reports the original data and (B) our updated values, if any. 
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'"'The corrected data used updated nonlinearities. ' fc 'The corrected data adjusted for cascading contributions. ' c 'The 
parenthesis indicates that the angle was not reported, so the value shown was the angle we believe was used. I"" The 
repetition rate was lowered to 760 KHz with the same result. MRt using quadratic nonlinearities of [fO]. "'The THG 
tensor component was measured instead of the cubic self-action components. ' g 'Fit using quadratic nonlinearities 
of ||||. W The measured mixture was |cn cos 2 8 +cij sin 2 0. 



(p9[), all suitably converted to 850 nm with Miller's rule, we can from the corrected value of 
= 4.08 ± 0.57 • 10~ 22 m 2 /V 2 and Eq. @ calculate 



C33 



13.8 ± 15.8 -10" 22 m 2 /V 2 , A =0.85 11m (31) 



An exciting consequence of this result is that it points towards a negative, self-defocusing Kerr 
nonlinearity. However, as indicated the particular value is very uncertain (the uncertainty was 
calculated using the error propagation rules), which mainly stems from the low prefactor in 



front of the C33 term in Eq. (30). A part from that, it also relies on three other separate measure- 
ments, making it very sensible to the en, cio and ci6 values used. Here we mention that the 
large uncertainty of cio reported in fll^ ] plays a role; with the various possibilities reported 
there for the cio value we always get a negative C33, and its value can vary by a factor of two 
from that reported in Eq. (|3lj). 

4.4. Summary of experiments 

A summary of the data reported in the literature along with our corrected or updated values 
is shown in Table [I]. This table also gives an overview of the crystal angles, pulse duration, 
repetition rates, wavelengths and the tensor components accessed in the measurements. 

Our analysis of the experiments [|l4|-^o|] points towards using the Kerr nonlinearities sum- 
marized in Table ||. Only the susceptibilities are reported in order to minimize the errors when 
using the values for composite nonlinear indices. We have there also indicated the Miller's delta 



Table 2. Proposed nonlinear susceptibilities for the BBO anisotropic Kerr nonlinearity. 
Note the value is deduced from the C] i and c'io coefficients, and the C33 value is deduced 
from the c\ \ , cio and c\§ coefficients. The Miller's delta Ayj; are calculated from Eq. (poj). 





A 
[Mm] 


x (3) 

[10- 22 m 2 /V 2 ] 


[10- 24 m 2 /V 2 ] 


Ref. 


(3) 

c n — Xxxxx 

(3) 

c io — Xxxyz 
c 16 = Xxxzz 

(3) 

c 33 — ZzZZZ 


1.053 
1.053 
0.850 


-0.24 ±0.04 
1.46 ±0.34 
-10.4±14.3 


52.8 ± 8. 4W 
-3.04 
23.4 
-286 


Eq.© 



'"'This value corresponds to en = 5.03 ■ 10 22 m 2 /V 2 at 800 nm. '*'The THG tensor component was measured 
instead of the cubic self-action components. 



calculated from Eq. (|26|). Since the experiments are performed at different wavelengths, using 
the AjjM values makes it easier to evaluate a linear combination of the nonlinear coefficient at 
some particular wavelength. [Another popular model for scaling the cubic nonlinearity is the 
Boling-Glass-Owyoung (BGO) model [|48|], but we did not see a big difference in using that 
model compared to Miller's rule. Besides, extending the BGO model to anisotropic nonlinear- 
ities is not straightforward, so we prefer to use Miller's rule for frequency scaling.] The values 
are reported with more significant digits than supported by the uncertainties, but this is done on 
purpose so it is easier to cross-check the coefficients as well as frequency scaling them. 

4.5. Implications for cascaded pulse compression in BBO 

In context of cascaded quadratic nonlinearities by far most important component is the FW SPM 

(3) 

coefficient j eff (©1 ; 0)i ), and in a type I 00 — >• e configuration it is given by the c\ \ component. 
Instead in a type ee — > e interaction it is given by the expression from Eq. (|Uj). In BBO and 
other borates this interaction is not so relevant as the quadratic nonlinearities are too small to 
support it, but in niobates it is extremely relevant as they are known for a large ^33 component: 
in this case d e s is maximized with = k/2, making C33 the effective cubic nonlinearity. 

In some of our previous work related to BBO ||3~7|, p9[|50|], we assumed an isotropic Kerr 
nonlinearity, and have used the value n 2 Kerr (G>i ; (£)\ ) = 3.65 • 10~ 20 m 2 /W taken from [p7|], and 
in a later publication, we have used the value given by Eq. (j25|). An implication of this larger 
value is that the "compression window" J3~7| ] becomes smaller. With this we imply a range of 
phase-mismatch values, where compression is optimal. We first need Ak to be small enough 
so «2 tot < 0, i.e. \n 2 casc | > n 2 Ken ((Oi; ©1). This is the upper end of the window, also denoted 
the Kerr limit. Once below the Kerr limit, decreasing Ak strengthens the total defocusing non- 
linearity. However, at a certain point the group-velocity mismatch (GVM) becomes too strong 
and the compression quality is strongly reduced. Several things happen: (a) the GVM-induced 
self-steepening term [|],|5j]] is increased as it scales as dyij Ak, where dyi is the GVM param- 
eter, so the compressed pulse experiences a strong pulse-front shock, (b) When Ak becomes 
too low (specified more accurately below) the SH will experience a resonant phase matching 
condition of a sideband frequency to the center frequency of the FW spectrum []52[], This is 
damaging to pulse compression because in the cascading process effectively this is the cas- 
cading bandwidth felt by the FW, and in the transition from the nonresonant to the resonant 
regime, the bandwidth essentially goes from being octave spanning to becoming resonant and 
thereby narrow [f26|]. The criterion for being in the nonresonant regime, also denoted the sta- 
tionary regime, can in the simple case where only second-order dispersion is considered [ |49[ ] 
be expressed as d\ 2 — Zk^ (oy^Ak < 0, where k^ (©2) is the group-velocity dispersion (GVD) 




Fig. 3. (a) Compression diagram for BBO type I cascaded SHG. In order to excite solitons 
Ak must be kept below the Kerr limit (red line). Optimal compression occurs when the 
cascaded nonlinearities dominate over GVM effects (Ak > M sr , above the black line). Note 
that Ak sr is calculated for the full dispersion case, and that for X\ > 1 .49 flm the FW GVD 
becomes anomalous. We have also indicated the operation wavelengths of Cnforsterite, 
Yb and TkSapphire based amplifiers. The Kerr limit employs Miller's rule to calculate the 
nonlinear quadratic and cubic susceptibilities at other wavelengths, case (1) corresponds 
to the 'old' Kerr value n? 2 Ken .(fi)i; £»i) = 3.65 • 1(T 20 m 2 /W@850 nm (taken from @) 
and case (2) corresponds to the Kerr value proposed in this work Ami = 52.8 rrr/V 2 , 
corresponding to KelT ((0i ; fi>i ) =5.1- 10~ 20 m 2 /W@850 nm. (b) Numerical simulation 
of the case marked with 'x' in (a): a 50 fs@ 1030 nm L m = 500 GW/cm 2 pulse propagating 
in a 25 mm BBO crystal with d = 19.1° and <j) = -90° (Ak = 80 mm -1 ). The simulations 
used the SVEA equations (^)-(^7|) including full dispersion and extended to include self- 
steepening, and case (1) assumes an isotropic Kerr nonlinearity and n\ KelT (ft>i \(0\) = 3.65 • 
10~ 20 m 2 /W@850 nm, while (2) uses the anisotropic coefficients of Table 



coefficient of the SH. Thus, in the case where the SH GVD is normal (positive) we have the 
stationary (nonresonant) regime with broadband cascading when Ak > Ak sr = d\il\Lk\ (o>i)]. 
Below this threshold, in the so-called nonstationary regime, the cascaded nonlinearity is as 
mentioned resonant: the poor bandwidth implies that there is no possibility to achieve few-cycle 
duration [[30|]. At the same time the compressed pulse quality is low as there is a strong pulse 
shock front: this stems from the dnjAk ratio being large, and this term controls the cascading- 
induced self-steepening, as mentioned above. As a consequence one can only achieve a decent 
pulse quality by keeping very low soliton orders (roughly below 2) [[30||. The implications of 



the larger nonlinear SPM Kerr nonlinearity for the FW is visualized in Fig. ^|(a). With the 
lower value, «2Kerr( ft) i ;fi) i) = 3-65 • 10~ 20 m 2 /W@850 nm, a large compression window is 
predicted, almost even encompassing TkSapphire laser wavelengths. With the new larger value, 
Ann = 52.8 m 2 /V 2 corresponding to n ! 2 KeTr {(0i; COi) = 5.1 ■ 10~ 20 m 2 /W@850 nm, the com- 
pression window starts opening around 900 nm. To judge the impact of using this larger Kerr 
value, Fig. ||(b) shows a simulation at 1030 nm, where the phase mismatch is chosen to lie 
inside the compression window using the new larger Kerr value. After 25 mm of propagation 
the short (50 fs) and intense (/;„ = 500 GW/cm 2 ) input pulse is compressed somewhat due to 
the formation of a self-defocusing soliton, and the spectrum shows SPM-like broadening and 
even formation of a soliton-induced Cherenkov wave (dispersive wave) around 2900 nm [JI3]. 
When instead using the 'old' lower Kerr value (dashed lines) the result is changed quite a lot: 
this is because |«2,tot| is now much larger so the intensity gives a larger effective soliton order 
Af 2 ff « |«2.tot|^in (the specific numbers are N e s = 3.6 in case (1) and A^ff = 1.9 in case (2). The 
soliton is therefore compressed more and earlier in case (1), so after 25 mm soliton splitting 
has already occurred and in the spectrum the broadening is more pronounced and the dispersive 
wave is much stronger. These results clearly show how sensible cascading is to the value of the 
FW Kerr SPM parameter. 

The experiments carried out at 800 nm were performed well into the nonstationary 
region, which was necessary to achieve a defocusing nonlinearity. No few-cycle compressed 
solitons were observed in [R], as the intensity had to be kept low in order to avoid severe GVM- 
induced self-steepening. (Note that the simulations in |Q] used = 5 • 10~ 22 m 2 /V 2 @800 
nm, i.e. identical to the value we suggest.) Instead the experiment carried out at 1250 nm [ |To[ ] 
was done in the stationary regime, and indeed a few-cycle soliton was observed. 

A question is whether pumping with e-polarized light could give any advantages. As the 
ee — > o interaction is heavily phase mismatched in a negative uniaxial crystal, the main cas- 
cading channel would be the ee — >• e interaction. As expected this interaction can never be 
phase matched. The maximum d e f[ is for 0=0, and this turns out to give the maximum cas- 
cading strength as well. However, our calculations show that it is too small to overcome the 
n 2 Kerr( ft) i ; ®i) nonlinearity, and the total effective cubic nonlinearity would be focusing. If 
9 = n/2 we get d e ff = da, which for BBO is very small (this is typical for borates, while the 
niobates instead have very large c/33 components). Nevertheless, for 9 = n/2 the C33 compo- 
nent determines the SPM coefficient, and as we found it could be negative. The total effective 
cubic nonlinearity would therefore be de-focusing. Investigating a BBO crystal pumped with 
9 = n/2 using e-polarized light would therefore be an interesting next step to resolve this issue. 

We should finally mention that the error made in assuming an isotropic response for the 
cascaded soliton compression is not large: most importantly, the crucial parameter is the FW 
SPM coefficient, which as we saw is the same in the isotropic and in the anisotropic cases for 
this type I interaction. The XPM and SH SPM coefficients change when anisotropy is taken into 
account, but they play a minor role and rather tend to perturb the result more than shape it. In 
order to make a more quantitative statement, we calculated the 9 -ranges where self-defocusing 
solitons can be excited in BBO. A necessary (but not sufficient [|37[]) criterion for this to happen 
is Ak > as well as AA: low enough for \n ! 2 casc | > n 2 Kerr The ran 8 e i s shown vs. FW 

wavelength in Fig. [|(a). Inside this 9 range we can now for each wavelength calculate the 
anisotropic Kerr nonlinearities. As we know the FW SPM coefficient does not change with 9 
(but it does change across the wavelength range shown, which we here estimate using Miller's 
rule), so we can normalize the results to this value. We then get the results shown in Fig. [Jb). 
The XPM term is seen to lie close to the isotropic value 1/3 (found by taking %jQ (a>i ; CO2) = 
en/3, dashed grey line). Instead the SH SPM term lies well below the isotropic value [found 
taking % eS ' ((O2; 02) = en, dashed blue line]. The main reason for this is the large negative C33 




\ [nm] l t [nm] 

Fig. 4. Estimating the operating parameters vs. FW wavelength for cascaded SHG in BBO. 
(a) The 6 -range for which Ak > is achieved as well as | casc | > n ! 2 Ken . ( a>i ; ft>i ) . (b) The 

XPM term X&[(®1 '• C0(2) and SH SPM term ^ (tt)2; fih), both normalized to the FW SPM 

term %jg {(0\ ', Oh). The range indicated corresponds to the 6 range in (a). The dashed lines 
indicate the isotropic limit. All nonlinear coefficients are scaled to other wavelengths using 
Miller's rule, and we used the SHG nonlinearities at 1064 nm from App. ^| as well as the 
cubic nonlinear parameters listed in Sec. a. We also took <j) = — 7C/2. 



component. In summary, assuming an isotropic response in BBO is a good approximation for 
the FW XPM term, mainly because the working 6 range lies close to zero (the isotropic limit). 
This result cannot be generalized to other nonlinear crystals, where the working range might lie 
closer to 7T/2; an example is lithium niobate that operates quite close to 9 = n/2 in the type I 
configuration, see e.g. [3^]. For the SH SPM term its value is quite far from the isotropic case, 
but off all the parameters this is the least important one because the SH SPM is negligible in the 
cascading limit compared to FW XPM and in particular FW SPM. To confirm this we checked 
that the simulations of case (2) in Fig. ||(b) were almost identical results when using isotropic 
Kerr nonlinearities. A similar conclusion was drawn in our recent work, where the anisotropy 
was also taken into account [^TJ . 

5. Summary 

We have analyzed the Kerr nonlinear index in BBO from a number of experiments [[l4|-|2"o[], 
and we argued that in many of them contributions from cascaded SHG nonlinearities need 
to be taken into account. The corrected nonlinear coefficients point towards using the Kerr 
nonlinearities summarized in Table |. They encompass all four relevant tensor components for 
describing the anisotropic nature of the BBO Kerr nonlinearity under Kleinman symmetry. The 
most reliable value is the cn component, relevant for describing the 0000 SPM coefficient, 
while the other three are more uncertain: the cig and cjo values are measured for THG, so using 
them for describing SPM and XPM effects is an approximation. Finally, the C33 coefficient is 
very uncertain as it was measured in a non-ideal way (since the experiment used a crystal angle 
where its importance is very small) and on top of that all the three other coefficients in the 
table were used to deduce its value. In all cases the measured coefficients had contributions 
from both instantaneous electronic and delayed Raman effects, but we assumed that the Raman 
contributions to the measured nonlinearities were vanishing. 

The proposed c\\ value was calculated as a weighted average over 10 different experiments 
in the 532-1064 nm range, and is therefore expressed as a wavelength-independent Miller's 
delta Ami = 52.8 ±8.4 x 10~ 24 m 2 /V 2 , cf. Eq. (§7|); at 800 nm this value corresponds to 
cn =5.03 x 10~ 22 m 2 /V 2 and n\ Kerr = 5.15 • 10~ 20 m 2 /W, which is larger than typical values 
used in the literature. When plotting the (corrected) values vs. wavelength, see Fig. g, we came 



to a surprisingly good absolute agreement with the 2-band model [[23[], except in the short end 
of the visible range. 

Finally we showed the predicted consequences for cascaded femtosecond pulse compression 
exploiting type I interaction in BBO when using this larger value: the operation range where 
few-cycle pulse compression can be achieved will shift to longer wavelengths, roughly above 
1.0 /im. We also showed that the anisotropic XPM coefficient of the FW is quite similar to the 
isotropic one, but this relies on the particular range of phase-mismatch rotation angles used in 
BBO, so this result cannot be generalized to other crystals exploiting type I interaction. Instead 
the anisotropic SH SPM coefficient is quite different from the one found assuming an isotropic 
Kerr nonlinearity, but since in cascading the SH is weak this difference will amount to much in 
a simulation of cascaded SHG. 
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A. The propagation equations under the slowly-varying envelope approximation 

We have previously touched upon the issue of anisotropic cubic nonlinearities in quadratic 
nonlinear crystals [|3(]], where we focused on type I interaction and the measurements in the 
literature for the lithium niobate crystal. Here we formulate the complete propagation equa- 
tions that hold for both type and type I, and discuss the anisotropic cubic tensor components 
relevant for the BBO crystal class. 

Let us show briefly how to derive the basic slowly-varying envelope equations (SVEA) for 
degenerate SHG (i.e. where the FW photons are degenerate), thus describing both type and 
type I SHG. We do this in order to avoid confusion about how to define the relation between the 
electric field and the susceptibility tensors. Working in mks units, from Maxwell's equations 
and applying the paraxial approximation we immediately get the time-domain wave equation 

V 2 E(r,f) = ^ ^2D(r,f) (32) 

where D = £oE + P is the displacement field, and P is the induced polarization. In frequency 
domain, where E(r, co) = /_!° 00 dfe" a 'E(r,f) and similarly for P(ffl), we may separate the linear 
and nonlinear polarization response as 

P(r,t») = £oZ (1) M -E(r,ffl) +P NL (r,co), (33) 

where ^'(ft)) is the linear susceptibility tensor. The nonlinear part is usually expanded in a 
power series as 

e^P^r, co) = x {2) (co) : E(r, e>)E(r, co) + x} 3) (fl»):E(r, ffl)E(r, o)E(r, co) + ■ ■ • (34) 

The tensor products are calculated by summation over the indices £q 1 Pnlj = HjkXijk^j^k + 
Hjkl XijkiEjEkEi H .By assuming plane waves (no diffraction), we have V 2 = Jp- so 

d 2 CO 2 CO 2 

^E(z, co) + ^e(co) -E( z , co) = -_ef 1 p NL ( z , co) (35) 



where £.(©) = 1 + % C®) i s me relative permittivity tensor, and we used £o/io = 1 /c 2 . 

For SHG we define the electric field envelopes Sj for the FW (frequency C0i) and the SH 
(frequency <x>2 = 2a>i), as well as envelopes for the nonlinear polarization response ^nlj 
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ui < ri(z,0e^ l(£ai)z " , ' £ai ' + u 2 ^ 2 (z,0e , ' ,:2(fi,2) ^" i>2 ' + < 



Pnl(z,0 = t [u 1 ^ NLl i(z,?)^^) z - tol '+u 2 ^ NL)2 ( Z ,0^ 2( ^ )z ^ + c.c. 



(36) 



(37) 



where the wave numbers are defined as k 2 (co) = p-£/((») with £j(co) = Uy •^W(fi)), and where 
u is a unit polarization vector. In a lossless medium £; is real, and the linear refractive index is 
given by rij((o) = weJ(co). The envelope wave equations are therefore 



dz 2 
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[kj(a>)-kj(a>j)]£j(z,(a) 
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The SVEA assumes |^| < kj(a>j)\gj\ giving 



d<gj{z,(0) r , . . 

1 ] y z J +[kj(«>) 



■kj(aj)]£j(z,a>) 



CO, 



2nj((Oj)c 



e ^NLjfcCO) 



(39) 



where we used the approximation k 2 (co) — k 2 ((Oj) — 2kj( (Oj ) [kj ( ffl) — fc; ( CO/ ) ] . At this time, we 
have also used that the assumption of a slow envelope compared to the fast oscillating carrier 
term e _,a V means that the nonlinear polarization term on the right-hand side can be approx- 
imated as f 1 , .ij Co^-^NL /(©) — . s .j £rT 1 ^nl. /(to). We can now perform the Taylor 
expansion of the wave number as kj(co) = ^ =0 /n! _1 fey (fflbO(fi) — fl)y), where the dispersion 
coefficients are defined as kj(co) = d" 1 kj((o)/d(o"' . Following this, we multiply both sides 
with e~' n ' where the frequency detuning is £2 = co — co 7 -, and perform an inverse Fourier trans- 
form over £2 to time domain. This gives 



CO; 



2nj(cOj)c 



e 1 &>nl, j(z,t) 



(40) 



This is the SVEA equations in the stationary laboratory frame, describing a pulse envelope 



propagating with group velocity 1 /kj(c0j). The time-domain dispersion operator is 



(41) 



and it can be truncated at any order to study effects higher-order dispersion. 

Besides a trivial transformation to a suitably chosen co-moving frame, what is left now is 
to evaluate the tensorial contributions to the nonlinear coefficients. Due to the vast variety of 
tensor combination possibilities for anisotropic nonlinearities it is here convenient to specify 
the problem very specifically: we therefore constrict ourselves to SHG in a uniaxial crystal, 
i.e. where light can be o-polarized with unit vector e" or e-polarized with unit vector e e [ |30[ 

Eq. (B6)], and we also restrict ourselves to degenerate SHG, where the two FW photons are 

(2) 

indistinguishable. Let us define the quadratic nonlinear polarization response as P N ^(z 7 t) = 
yy) ■ E(r,f)E(r,f ), where we have assumed that the quadratic nonlinearity is instantaneous in 
time. The procedure is then to insert the field envelope in this nonlinear polarization response 



and then calculate the polarization response for a specific field by applying the dot product with 



the fields unit vector: e l uj ■ P^(z,t ). The only relevant terms are those that oscillate with the 



,(2) 



fast carrier frequency oscillations of the field e ,m i', 
3? . Doing this we get for the FW and SH fields 



as this defines the polarization envelope 



Mkz 



,-iAkz 



(42) 



where the phase mismatch is Ak = £2(02) — 2k\ (ffli ). The effective susceptibility is calculated 
from the combination = uj-X^ ■ u^u/ that appears from isolating the e~ m i' term. The X^ 
is a rank 3 tensor with 27 different elements, but due to permutation symmetry in the the crystal 
axes indices there are only 18 nonzero elements [^3j. Depending on the crystal symmetry class 
(point group) many of these are zero. It is therefore fruitful to consider a specific point group, 
and perform the calculations. We show in Sec. [| the results for the point group 3m. We note 
that the nonlinear polarizations would have been a factor of 2 larger had we not used a factor 5 
in front of our envelope definition. 

Moving to the cubic nonlinear tensor component, it is done in much the same way, and we 
refer to [|o[ App. B] for details. We again consider an instantaneous response, and the cubic 
nonlinear polarization response that oscillates with e~ m i t turns out to be related to SPM and 
XPM; we could have had contributions from third-harmonic generation here if we had included 
such a harmonic field Si(z,t)e ik3 ^ z ~ i(0}t in the envelope definition. We then get [|(| Eq. (B8)] 
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c ^NL, 



(3), 



(43) 



(3) 

The notation of the cubic susceptibilities x^n ty) implicitly assumes phase-matched SPM 
and XPM interaction between 0), and C0j, as opposed to general four-wave mixing. We always 

have Xgff (©1; ©2) = X^s (®2> ®i )■ The factor ^ comes from the factor j in front of our envelope 
definitions, Eqs. @-@); if we do not have this factor 5, as in e.g. Boyd [^8|], it is important to 
note that the values (and consequently also the n\ Kerr -values) in the two cases remain the 

same. The x is a rar, k 4 tensor with 81 different elements, but due to permutation symmetry 
in the the crystal axes indices there are only 30 nonzero elements p3[]. The calculations of the 
Xtfl specific to a uniaxial crystal in the 3m point group was done in [ p0| App. B], and in Sec. [| 
we summarize the results. 

We can now express the plane-wave SVEA equations for the electric field envelopes as 
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(44) 



(45) 



The time x follows the FW group-velocity 1 /k^\(Oi) by the transformation from the lab time 
t as x = t — zk^\(0\ ), which gives the group-velocity mismatch (GVM) term dn = k^\c0i) — 



k^'(a>2)- We use the short-hand notation rij = nj(G)j), and af e ff = XeS ' s me usua l reduced 
notation of the effective quadratic nonlinearity. We can convert the electric field to intensity 
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Sj — > (2/eon J c) 1 / 2 Aj, so |A ; -| 2 is the intensity in [W/m 2 ], and we get 
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where the equations now have a common SHG nonlinear parameter 

/ _ (Q\d e « / 2 



(46) 
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We can also establish the link between the Kerr nonlinear refractive indices and the cubic non- 
linear susceptibilities as 



(49) 



cf. also Eq. (C6) in [|30|] . Both «2 Kea- an d z' 3 ' are in mks (SI) units, while conversion to and 
from the esu system is reported in Q30|]. 

The SVEA equations can be extended to the slowly-evolving wave equations (SEWA) |l 
by including self-steepening effects, and also the non-instantaneous (delayed) Raman effect can 
straightforwardly be included [p7|, pip . For simplicity such effects are neglected here. 



B. Reduced nonlinear Schrodinger equation in the strong cascading limit 

The cascading nonlinear contribution in the strong cascading limit (AkL 3> 2n) can quickly 
be found by taking assuming an undepleted FW. Thus E\ is taken independent on z in Eq. 
(p3]), and the SH can therefore be solved directly when the Kerr SPM and XPM terms are 
neglected (which follows directly from the undepleted FW assumption). In absence of disper- 
sion (the effect of dispersion is discussed in Sec. i.5) we get that the SH is slaved to the FW: 
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Sf{z, l)e ' Akz . When inserting this in the FW equation (44) we 
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arrive at the following nonlinear Schrodinger equation for the FW 
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(50) 



From these equations we now understand the basis of the nonlinear index change described by 
Eq. (0): in the strong cascading limit the FW experiences a cubic Kerr self-action nonlinearity 
| <5ip, where both cascading and material Kerr nonlinearities contribute. The cascading Kerr- 
like nonlinearity is described by the parameters 
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(51) 



where the latter is identical to Eq. (|lj) and it comes from using Eq. (|49|), or simply performing 
the same exercise with the intensity version of the SVEA equations. These will be used in Sec. 
to estimate the cascading contribution in the various experiments. A practical note: this lowest 



order cascading nonlinearity seemingly diverges at Ak = 0, but this is because the assumption of 
an undepleted FW breaks down in this limit. To the next order a correction term 1 — sinc(AfcL) 
must be applied [^3j, which resolves the divergence at Ak = and from this the maximum 
cascading strength is obtained at AkL ~ n. 

There is also a quintic nonlinearity from cascading mixed with the XPM term of the FW 

(5) _ ncofxiff\coua>2)d 2 en , 4(B 1 2 »{ tor (0)i;c^)<4 

* CaSC ~ 5«2 C 2 A £2 ' "4,casc- ^2^3 ^) 

where the latter comes from the relationship n\ = 5^^ 5 - > /(4«i£qC 2 ), and the quintic nonlinear 
refractive index n\ was defined in Eq. (13). This term gives a n online ar index change that is 



positive and °< I , and becomes important for high pulse fiuences [ |37| , p2fl . However, we evaluate 
it to be negligible in the experiments treated in Sec. ^. 

Note that these simple cascading expressions only hold in the strong cascading limit, which 
essentially means that the characteristic length for the phase mismatch (i.e. the coherence 
length) is much shorter than any other characteristic length scale in the system [^]. This can 
be achieved by fulfillment of two criteria: (a) a strong phase mismatch AkL ^ 2n (i.e. the many 
up- and down-conversion steps must occur inside the interaction length), (b) A phase mismatch 
that dominates over the quadratic nonlinearity. To express this we introduce the traditional [^] 
quadratic nonlinear strength parameter T = »ic/ e ff Kin|/(v /w i n 2c), where |<S; n | is the peak input 
electric field. Criterion (b) can then be expressed as Ak 3> T (which essentially ensures an un- 
depleted FW, see also the discussion in [|3]]). In all the cases in Sec. ^ where we correct for 
cascading contributions, we evaluate that the strong cascading limit is fulfilled. 



C. BBO crystal parameters 

The crystal beta-barium-borate (/3-BaB2C>4, BBO) is a negative uniaxial crystal of the point 
group 3m, which has a transmission range from 189-3500 nm [|29"|, |54|] . The Sellmeier equa- 
tions were taken from Zhang et al. [|43|]. We used the quadratic nonlinear tensor components 
measured by Shoji et al. [||]: at A = 1064 nm \d 2 2\ = 2.2 pm/V, |d 3 i| = |<f 22 |0.018 = 0.04 
pm/V, d l5 = 0.03 pm/V, and d i3 = 0.04 pm/V; at A = 852 nm \d 22 \ = 2.3 pm/V; at A = 532 
nm \dzi\ = 2.6 pm/V; at A = 1313 nm \dzz\ = 1.9 pm/V. Note that Kleinman symmetry has 
<^3i = d\s, which is only slightly violated at 1064 nm by these measurements; we therefore 
throughout this paper assume Kleinman symmetry. The sign of the product d^\d 22 has been 
shown to be negative J55|], which means that flipping the crystal 180° gives a changed effec- 
tive nonlinearity (see also discussion for Fig. |]). An overview of other measurements is given 
in J56|]; in particular note that the 1064 nm d 22 value is similar in most measurements, but the 
d\s value was historically measured to be higher. This stems from the initial measurements 
of BBO [|4]] that gave d\ 5 = Q.01d 22 [§|], which then with the d 22 = 2.2 pm/V value gives 
d\5 ~ 0.15 pm/V. A users note of the Shoji et al. values: the 852 and 1064 nm measurements of 
the \d 22 \ value agree well with Miller's scaling (they have the same Miller's delta |p2[|), while 
the 1313 nm measurement seems to lie too low. Therefore when calculating the effective non- 
linearity at an arbitrary wavelength the safest option generally is to use the 1064 nm values, 
where all the nonlinear components were measured, and employ Miller's scaling to go to other 
wavelengths. 



